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We examine the transmission through nonideal microwave resonant circuits. The 
general analytical resonance line shape is derived for both inductive and capacitive 

q_; 

coupling with mismatched input and output transmission impedances, and it is found 
that for certain non-ideal conditions the line shape is asymmetric. We describe an 
analysis method for extracting an accurate internal quality factor (Qi), the Diameter 
Correction Method (DCM), and compare it to the conventional method used for 
millikelvin resonator measurements, the <fi Rotation Method (</>RM). We analytically 
find that the 0RM deterministically overestimates Qi when the asymmetry of the 
resonance line shape is high, and that this error is eliminated with the DCM. A 
consistent discrepancy between the two methods is observed when they are used to 
analyze both simulations from a numerical linear solver and data from asymmetric 
coplanar superconducting thin-film resonators. 
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I. INTRODUCTION 



Precise measurements of both loaded and internal quality factors of thin-film supercon- 
ducting resonators are necessary for many applications, from astronomy photon detectors,- 
to materials analysis, 2 '- to qubit readout.-^ Growing interest in superconducting quantum 
computing has recently motivated detailed measurement of several types of superconducting 
thin-film resonators at millikelvin temperatures.-^ - — Unfortunately, non-ideal experimen- 
tal setups can lead to an asymmetry in the resonance line shape, -^^i^ - — corresponding 
to a rotation of the resonance circle, which complicates the interpretation of internal and 
external quality factors (Qi and Q e ). Several methods for extracting the Qi of a resonator 
exist for different experimental setups.— - — However, these methods either require a single 
port reflective measurement^^ (incompatible with most qubit measurements), full two-port 
data^ 1 ^ (typically unavailable for millikelvin measurements), or identifying and fitting to 
a second coupled modo^ (a special case). In contrast, the most widely used technique for 
analyzing millikelvin resonator measurements, the <fi rotation method (0RM), simply adds 
an empirical rotation of the resonance circle to extract the quality factor.-^^ 

In this article we show how asymmetry in the resonance line shape can arise from cou- 
pling the resonator to mismatched input and output transmission lines and non-negligible 
transmission line series inductance. Based on this understanding of the origin of the asym- 
metry, we derive the Diameter Correction Method (DCM), used in recent publications,— 1 ^ 
and use it to extract Qi. We compare this to the conventional analysis method, the 
rotation method (0RM)^^, and show that there is a one-to-one mapping between the 
two methods but that the 0RM systematically overestimates Qi by an analytically quan- 
tifiable amount. Note, we will not address fitting techniques here because a comprehensive 
quantitative comparison of fitting techniques has been made.~ 



II. DERIVATION OF ASYMMETRIC RESONANCE 

We consider a notch type resonator coupled to input and output transmission lines (see 
Fig. DJa)) in which the transmission, S 2 \ = V out /V in , is measured. Full transmission is 
measured off resonance and reduced transmission is measured on resonance. The resonator 
inductance and capacitance are L and C , where C is complex to account for dielectric 
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loss. Vi n and V^t are the input and output voltage waves, and V is the voltage across the 
capacitor, C. Ideally, the transmission line ports are matched (Z in = Z out = Z ), and L\ is 
small (L x << CcZf n and L\ « CcZ% ut ). However, conditions such as a precise matching 
of Z in and Z out can prove to be difficult to achieve in an experimental setup. 

Assuming a high quality factor, Qi » I, the circuit in Fig. [TYa) can be redrawn as 
Fig. [T](b), where R = Qi/ (ojqC). Solving Kirchhoff's equations we find an expression for the 
transmission as a function of the voltage across the capacitor: 

s - = (1+ K 1 + ^(t +z ^ Cc ))- (1) 

where 1 + e = — -, — ■ — ^ — , Z'- = Z in + iuL 1 — iu^-, and lei << 1. Solving Kirchhoff's 

equations again and eliminating V out we get the independent equation: 

V\ —+iuC+ , ° + 

zcjL tujCcZ out + 1 

A/ iuiCcZput 

L 1+iuiCc Zout 



Z' out + Z in + iu (Li-^) 



L l+iwCcZout 



ml ,^ t + ^ + ^(L 1 -Mi) y l' {2) 

where Z' out = - • Equation (JSJ) is of the form V ( ^ + iwC + Gn) = In, where Gn 

and In are the Norton equivalent conductance and current respectively. The circuit thus 
has a Norton equivalent shown in Fig. Hfc). From Eq. (j2j) we see that 

/ M iuiCcZout \ 

W = -^^m | — — — ; —77 ) > W 



and, 

M iuiC'cZ ou t 



iuCc \L 1+iwCcZouty f |, 

N = iuC c Z out + 1 + Z' out + Z m + zoo (Lx - 5) ' U 
Since the real part of G^r loads the resonator measurement and the imaginary part shifts 
the resonance frequency, it is useful to separate Gn into its real (conductive) and imaginary 
(susceptive) components. We define RelG^} = 1/Rt, Imi^G^} = ojCt, R~tt = R~ x +R^ , 
Q = R e ffU (C + Ct), and G' = (77 + Z' in iuCc)- Then, Eq. jT]) can be rewritten as 
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We now note that to second-order in the small parameters (M/L) and uCcZ out , G' is 
equal to . Stopping at second-order would yield a resonance with a symmetric Lorentzian 
line shape. To expand to higher-order, we rewrite Eq. (j^D as 



where 



G £) = G' — Rj,^ (7) 



Expanding to third-order we find: 



G D = iuCc^- ( Z J n " f "* I + 



,2 



(f)V-^L)), (8) 



and note that Go is purely imaginary. We now define 

w (C + C r 



and recognize that for Ct « C 

Qr^g-i-RejQ; 1 }, (10) 

Thus we find 

g 8 l = (l+g)fl- i ,^V, )■ (11) 



1 + 2zQ^ 



We emphasize here that Go is the term that creates the asymmetry in the line shape 
or equivalently a rotation of the resonance circle around the off resonance point. If Go 
were zero then Q~ l would be real, reducing Eq. ffTTj) to a symmetric Lorentzian line shape. 
Equation (TJTJ can be written as 



Q Qe 1 

S 21 = (l + e) I 1 



1 + 2iQ 



(12) 
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where Q e 1 is represented in terms of its magnitude and phase, <fi. Another equivalent 
representation is 



Qe \ k>0 

1 + 2iQ^^ 



S 2l = (l + e) , (13) 



where we have defined 1/Q e = He ^Q^ 1 ^ and 5u is the difference between the resonance 
frequency and the new rotated in-phase point on the resonance circle, U\ (see Fig. |2j). The 
form of Eq. f lT3"j) can be understood by noting that S21 is real when ou = Sou + Uq (to within 
a phase rotation of 1 + e). 

Here we stress that Eqs. ( fHlfTBl) are equivalent representations of the asymmetric line 
shape, each highlighting a different interpretation of the asymmetry. In Eq. (llip the asym- 
metry is quantified by Im jQ^j and one can think of the asymmetry as coming from a 
complex loading of the resonator. In Eq. (1121) the asymmetry is quantified by 0, where is 
the rotation angle of the resonance circle around the off- resonance point (see Fig. [2]). And 
finally in Eq. ffl3|) the asymmetry is quantified by Sou, where 5u = ou\ — Uq is the frequency 
shift of the in-phase point on the resonance circle from uuq to ou± (see Fig. [2]). There of course 
exists a simple one-to-one mapping between the three notations: 

/ Im{Q-n\ / Sou 



(b = arctan ^ — = arctan 2Q — . (14) 

\Re{Q 7 i}J V V ^ ; 

One method that has been used^ 1 ^'^ to extract internal quality factors simply accounts 



for the asymmetry by adding an empirical rotation and incorrectly substitutes 
Q~ x by defining 



for 



1 1 



1 



Qe 



(15) 



Qi, (j>RM Q 

This method accounts for the asymmetric line shape phenomenologically by adding the 
rotation, 0, without accounting for its origin and its impact on the interpretation of Qi. 
It corresponds to rotating the resonance circle back an angle 0, thereby putting ouq on the 
in-phase axis. This is the <fi rotation method (0RM), and the rotation of the 0RM can 
best be seen by examining the difference between Fig. [2^a) and Fig. [2(b) (A). However, 
simply rotating the resonance circle by angle <fi does not take into account the fact that the 
asymmetry has also caused the circle to grow by a factor of 1/ cos(0), assuming the circle 



has been normalized to full transmission off resonance (S , 2i(w << coq) = S21 » <^o) = !)• 
We have shown here that instead one has 

1 ~ 1 1 (16) 

Qi, DCM Q Qe 

We call this the Diameter Correction Method (DCM) because in addition to rotating the 
circle by the asymmetry angle, 0, it also corrects the diameter by accounting for the com- 
plex Q e . This can be seen by examining the difference between Fig. [2]^a) and Fig. [2]^b)(n). 
Another interpretation of this result is that the quantity that remains constant in the asym- 
metry transformation is not the diameter of the resonance circle, as the 0RM assumes, but 
rather the distance between the in-phase axis intercepts, shown in bold in Fig. [3 That 
invariant length is the diameter of the circle for a symmetric resonance and becomes a chord 
of the circle as asymmetry is added, but remains equal to Q/Q e , while the diameter grows 
as Q/\Q e \- The analytical discrepancy between the two methods can simply be determined 
by subtracting Eq. ( fl6l) from Eq. ( fl5l) . 



1 1 



1 



Qe 



(cos(0) - 1) . (17) 



Qi, DCM Qi, 4>RM 

From Eq. ( IT71) we see that the error in the 0RM diverges for high asymmetry angle, ~ ±7r, 
and for low Q e , high coupling. Note that for = 0, Q~ x is real and Eq. (Tl6|) reduces to 
Eq. ffT5"]) and therefore Eq. ffTTj) goes to zero. 



III. FITTING AND ANALYZING SIMULATIONS AND DATA 

To test the 0RM and the DCM, we simulate the transmission using a numerical linear 
solver. Simulations are run varying several different parameters: Qj, impedance mismatches, 
strength of both inductive and capacitive coupling and inductance of L\. The resonator ca- 
pacitance and inductance are held at 0.3 pF and 2.5 nH, respectively, producing a resonance 
frequency that ranges from 5.717-5.802 GHz (resonance frequency varies with coupling ca- 
pacitance). The simulated data is then fit and analyzed using both methods. 

We created asymmetry by varying Z in /Z out . Note that asymmetry can also be created 
by increasing L\. However, L\ values in the nanohenries are required to create significant 
asymmetry which is far too large to be physical. Figure [3] shows results from simulations 
and fits with the same simulated quality factor, Qi = 10 5 and a range of Z in /Z out values. 
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The coupling line mismatch creates a clear asymmetry in the line shape which is quantified 
with the extracted asymmetry angle, <f), also shown in Fig. [3j In addition to the value of </> 
extracted from the fit, we also analytically determine the asymmetry using 



<P = arctan (^^) • (18) 

Two internal quality factors are extracted from these fits, one using the 0RM and the 
other using the DCM. In Fig. H] both extracted quality factors as well as the fit extracted 
asymmetry angle <fi are plotted against the predicted (f) from Eq. ( I18p . It is clear from Fig. H] 
that the DCM is more accurate as the asymmetry, 0, increases, and that the two methods 
agree for small asymmetry. 

We also compared both analysis techniques when asymmetry is held constant but Qi is 
varied, which we will show later corresponds to some experimental data sets. Figure |5] shows 
the fit extracted Qi from both analysis techniques as the actual simulation Qi is increased 
for two sets of simulations, one with low and one with high asymmetry. For low asymmetry 
(matched ports) both analysis techniques yield the correct Qi within the expected first-order 
error (Cq/C). However, for high asymmetry (mismatched ports), the 0RM yields quality 
factors that are systematically too high. For sufficiently high QiS, the 0RM yields negative 
QiS (this is why the asymmetric data analyzed by the 0RM appears to stop at large Qi in 
Fig. [5]). These unphysical, negative, QiS can best be understood by examining the circle 
plots in Fig. [2j As discussed earlier when there is a large asymmetry, in addition to being 
rotated, the resonance circle grows by a factor of 1/ cos(<ft) (assuming full transmission off 
resonance). Since the 0RM only rotates the circle back, it does not account for the increase 
in size, shown in Fig. Mjoi). So if Q ~ Q e (Qi » Q e ), the circle diameter is larger than 
1, almost crossing the y-axis. Rotating the circle using the 0RM causes the circle to cross 
the y-axis and this yields a negative Qi. In Fig. EJ^b) the 0RM analyzed simulation almost 
crosses the origin. This corresponds to the Qi = 8 x 10 5 simulation in Fig. |5) it is an example 
of a simulation data set with a Qi and asymmetry not large enough to create a negative Qi 
but still large enough to create a considerable discrepancy between Q it ^rm and Q^ dcm- 

To further evaluate both methods, we also analyzed data from a 5.75GHz coplanar alu- 
minum resonator on sapphire, measured at 30 niK in a dilution refrigerator. Figure [5(a) 



shows a picture of this resonator and a more detailed description can be found in Ref. 
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The 



resonator is measured by being mounted in a copper sample box and electrical connections 



are made with the Coplanar waveguide (CPW) using aluminum wire bonds. Figure E(b) 
shows an example of the measured resonance line shape and its fit for one mounting of the 
resonator which exhibited a particularly high asymmetry, presumably due to the non-ideal 
mounting of the device in the sample box. In Fig. [7J we show the extracted Qi using both 
techniques and the extracted asymmetry angle, <fi. The quality factor dependence on voltage 
is discussed in Ref. [l0|. Here we focus on the difference between the two analysis techniques. 



Figure [7J is very similar to the high asymmetry simulations in Fig. |5j As expected from 
Eq. (|Sj), the asymmetry, 0, is independent of Qi for both the real device measurements and 
the simulated data. Also the last two data points for the </>RM in Fig. [7J are negative (and 
off the plot) in the same manner that the last points in the simulated data of Fig. are 
negative. 

An additional way to test the analysis techniques is by varying Q e while keeping Qi 
constant. In Fig. [HJ Q e is increased by lowering the capacitive coupling. As expected, 
for low asymmetry (matched ports) both analysis techniques do a good job of extracting 
Qi = 10 5 . However, with high asymmetry (mismatched ports) the 0RM overestimates Qi 
by a decreasing amount as Q e /Qi increases. Interestingly, in the 0RM, as Q e increases, 
the extracted Qi approaches the real value, although the asymmetry, 0, is increasing. This 
is because as Q e increases the </>RM is less sensitive to the asymmetry. This behavior is 
captured in Eq. (TTTj) . which shows that as Q e increases, the difference between the two 
analysis methods vanishes. In fact for Q e » Qi, the asymmetry becomes completely 
irrelevant and the two techniques converge. 



IV. CONCLUSION 

In summary, we derived an analytical resonance line shape based on circuit parameters 
and found that for non-ideal conditions the line shape is asymmetric. We developed a tech- 
nique (DCM) for extracting accurate internal quality factors from asymmetric resonator 
measurements using only transmission data. By analyzing simulated resonator measure- 
ments we found that the DCM is superior at extracting accurate internal quality factors to 
the conventional 0RM used in millikelvin resonator measurements. We found that in the 
limit where the asymmetry is low, the two methods agreed, but when the asymmetry is high, 
particularly when Qi » Q e , the DCM accurately determines Qi while the 0RM systemat- 
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ically overestimates it. Also, for sufficiently high asymmetry and coupling the 0RM gives 
a negative Qi. We have also shown that the two methods can produce different results on 
real data taken on a coplanar superconducting aluminum resonator with high asymmetry. 
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FIGURES 
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FIG. 1: (a) Schematic of resonator measurement setup with both inductive and capacitive 
coupling and mismatched transmission lines, (b) equivalent circuit to (a), where C has 
been separated into its capacitive part (C) and its resistive parts (1/uR). (c) Norton 
equivalent circuit for resonator measurement where V is the voltage across the capacitor, 
C, and G N = R T + l/uC T . 
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FIG. 2: (a) Simulated transmission through mismatched coupling lines (Z in = 24.5S1, 
Z out = 84.511) plotted as ImjS^i} vs. Re{S , 2 i} with a fit to a circle. The asymmetry is 
represented as a rotation of the resonance circle by the angle away from the real 
(in-phase) axis or equivalently as 5u = ui — u , the frequency shift of the in-phase point on 
the resonance circle, (b) Shows the simulated transmission with asymmetry removed using 
both analysis techniques. The 0RM (A) only rotates the circle to the real axis while the 
DCM (□) both rotates the circle and removes the factor of 1/ cos(0) increase to the 
diameter. The DCM shows that the invariant quantity is not, as the 0RM assumes, the 
diameter of the circle (equal to Q/\Q e \ and Q/Q e before and after the DCM 
transformation respectively) but rather the length of the real axis segment intersecting the 
circle (shown in bold and equal to Q/Q e ), where 1/Q e = Re 1 1/Q e > . 
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FIG. 3: Simulated and fit to symmetric and asymmetric resonance line shapes. Here the 
asymmetry is created using mismatched coupling lines (Z in and Z out ). Asymmetry angles, 
0, are extracted from the fits. 



14 



(3) HO. 

□ DCM 
. 0RM 



B- □ □ ft ti 




(j)pred. (radians) 



FIG. 4: (a) Qi extracted from fits to circuit simulations using both analysis techniques, 
0RM (•) and DCM (□), as a function of predicted asymmetry angle, (p pre d., calculated 
using Eq. ([TBI) . The dashed line indicates actual simulation Qi. At low asymmetry the two 
methods agree. As asymmetry is increased, the </>RM extracted Qi deviates from the actual 
Qi. (b) The fit extracted asymmetry angle, <p S i m , (♦), as a function of predicted asymmetry 
angle, (p pre d.- The solid line is the 4>pred. = 4> S i m . line. Good agreement of that line with the 
results (♦) indicates that this method is accurate at predicting the asymmetry. 
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FIG. 5: (a) Qi extracted with both analysis techniques (0RM and DCM) as a function of 
the actual Qi from two sets of simulations. The first set of simulations had high 
asymmetry (mismatched ports: Zi n = 24.5J1, Z out = 84.5f2) and the second had low 
asymmetry (matched ports: Z in = Z out = 50S1) . Solid line is actual Qi equal to extracted 
Qi line and dashed line indicates the coupling (Q e ). Both analysis techniques work well 
with low asymmetry but only the DCM works with high asymmetry at large Qi. At low 
simulation internal quality factors (Qi = 10 5 ) the DCM extracted internal quality factors 
(Qi = 1.002 x 10 5 ) with less than 1% deviation from the actual value in both low and high 
asymmetry simulations and at high simulation internal quality factors (Qi = 4 x 10 6 ) the 
DCM extracted internal quality factors (Qi = 3.85 x 10 6 ) with less than 4% deviation from 
the actual value for both low and high asymmetry simulations. The deviation at high 
internal quality factors is limited numerically by the fit and is not a limit on the method, 
(b) The fit extracted asymmetry angle, (p S im., for both low (0) and high (♦) asymmetry 
simulations plotted against the actual simulation Qi. 
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FIG. 6: An image of a coplanar superconducting aluminum resonator, coupled to a 
coplanar waveguide (CPW) transmission line. The resonator is composed of a meandering 
inductor (left) and a long coplanar strip (right). In one sample box mounting this device 
shows high asymmetry, (b) An example of a measured line shape of the resonator with its 
fit. Note that the fit is centered at the resonance frequency, but not at the minimum 
transmission frequency because those are not the same frequencies for asymmetric line 
shapes. 
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FIG. 7: Data from a resonator identical to that shown in Fig. 6. (a) Q iy extracted using 
both analysis techniques, 0RM (A) and DCM (■), as a function of voltage across the 
resonator, Vrms- As with simulated results in Fig. EJ the 0RM systematically extracts 
higher QjS with the highest QiS yielding negative results. The dashed line indicates the fit 
extracted external quality factor, Q e . (b) The fit extracted asymmetry angle, 4> exp . (♦), 
plotted against Vrms- Again similar to the simulated results in Fig. |5l (fi exp . is independent 
of the changing Qj. 



18 



(a) 

200 
^ 180 
1 160 

1 140 

£ 120 

<5 100 
80 
(b) 0.2 

? 
S, -0-2 
4 -0.4 
-0.6 



x 



□ 

A 
x 



matched (DCM) 
matched (0RM) 
mismatched (DCM) 
mismatched (0RM) 



X 



X 



_bA 



-bA- -sA bA- -bA 



















<0 

♦ 



matched 
mismatched 



10* 



10° 



Qe 



FIG. 8: (a) Qi extracted using both analysis techniques from two sets of simulations. One 
with high asymmetry (mismatched ports: Zi n = 24.5f2, Z out = 84.5f2) and one with low 
asymmetry (matched ports: Zi n = Z out = 50f2) plotted against a varying Q e . Q e is varied 
by varying the coupling capacitance (1-10 fF), with a constant mutual inductance (5 pH). 
The dashed line indicates the actual Qi of the simulations. With increasing Q e the 
inaccuracy of the 0RM is diluted due to the decreasing weight of Q e in the analysis, (b) 
The extracted asymmetry angles for the two simulations, low (0) and high (♦) asymmetry. 
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